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Abstract 

We test the Sine function representation, a novel method for numerically evaluating Feynman diagrams, by 
using it to evaluate the three-loop master diagrams. Analytical results have been obtained for all these diagrams, 
and we find excellent agreement between our calculations and the exact values. The Sine function representation 
converges rapidly, and it is straightforward to obtain accuracies of 1 part in 10 6 for these diagrams and with longer 
runs we found results better than 1 part in 10 12 . Finally, this paper extends the Sine function representation to 
diagrams containing massless propagators. 
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1 Introduction 

The Sine function representation, introduced by us in 
Ref. , is a new approach to numerically evaluating Feyn- 
man diagrams. This paper tests the Sine function rep- 
resentation by using it to calculate the three-loop mas- 
ter diagrams, whose analytical forms were obtained by 
Broadhurst in a mathematical tour de force . We have 
selected these diagrams because, in order to make a mean- 
ingful test of the Sine function representation, we must 
apply it to diagrams which contain multiple loops and 
which have been evaluated analytically, so that we can 
compare our results to known exact values. The criteria 
"topologically complex" and "evaluated analytically" are 
not mutually compatible, but the three-loop master dia- 
grams are sufficiently complex to allow a stringent test of 
our methods. In addition, the master diagrams contain 
massless lines, so we must start by extending the Sine 
function representation to theories with massless propa- 
gators. 

The Sine function representation performs an ab initio 
computation of the corresponding integral, so we make 
no use of the special analytic properties of the master 
diagrams. Consequently, our success with these diagrams 
supports our more general claim that the Sine function 
representation is a powerful tool for accurately evaluating 
arbitrary complex diagrams. 

We believe that our results, which are contain up to 12 
or 13 significant figures, represent an unprecedented level 
of precision for the numerical evaluation of a non-trivial 
Feynman integral. Conventional wisdom holds that multi- 
dimensional Feynman integrals must be numerically eval- 
uated with Monte Carlo methods, of which VEGAS is the 
canonical implementation || . Monte Carlo methods usu- 



ally yield the first few significant figures without too much 
difficulty but their efficiency drops extremely rapidly as 
the desired accuracy is increased, making highly accurate 
calculations effectively impossible. 

Obviously, in any realistic situation these diagrams 
correspond to higher order corrections, and we are capa- 
ble of a level of precision far beyond than that of any ex- 
perimental data. However, if only the first few significant 
figures of the exact result are need the Sine function rep- 
resentation can compute these very quickly. Furthermore, 
Broadhurst 's analytical treatment of the master diagrams 
relies on the PSLQ algorithm |4| for identifying integer 
relationships. Part of input required by this algorithm 
requires is accurate numerical evaluation of the diagrams. 
By first reducing the master diagrams to one- (or two-) 
dimensional integrals, Broadhurst was able to extract the 
numerical results he needed. However, if one wished to 
apply the same techniques to diagrams that cannot be 
reduced to numerical integrals which are tractable using 
standard methods, the Sine function representation pro- 
vides a promising alternative approach for their numerical 
evaluation. 

In the following section we summarize the Sine func- 
tion representation, and show how it can be applied to 
diagrams with both massive and massless propagators. 
In Section 3 we define the three-loop master diagrams 
and derive their Sine function representations. We evalu- 
ate these diagrams in Section 4 and compare our results 
with the analytic forms, while in Section 5 we discuss the 
implications of our work. 
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2 Sine Function Representation 

The Sine function representation |Q hinges upon a novel 
representation of the propagator and, as the name sug- 
gests, is derived using the properties of the generalized 
Sine function, 



Sk(h,x) 



sin [ir(x — kh)/h] 
tt(x — kh)/h 



(1) 



where the capital S distinguishes this from the more fa- 
miliar form, sinc(x) = sin(a;)/x. The properties of the 
Sk(h, x) are thoroughly discussed by Stenger ||. For our 
purposes, their most important attributes are that any 
function f(z) which is analytic on an infinite strip of the 
complex plane that includes the real axis can be approx- 
imated as 
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f(kh)S k (h,z) 



while 



f(z) dz 



k— — o 



f(kh). 



(2) 



(3) 



Crucially, the accuracy of the approximation improves ex- 
ponentially as h — ► [||. In Rcf. |IJ] we constructed the 
Sine function representation for the scalar field propaga- 
tor which, in Euclidean momentum space, has the form 



Ga(p) = 



1 



-P 2 M 2 



P 



(4) 



where we have introduced a cutoff in anticipation of di- 
vergent integrals. We now derive a version of the Sine 
function representation which is applicable to both mas- 
sive and massless fields. 

The spacetime propagator is the Fourier transform of 
G A (p), 



G A (x) 



1 



(4tt) 2 



ds 
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where we introduced s by exponentiating the denomina- 
tor, and then performed the Fourier transform. We now 
replace s by e z and use equation (|^) to approximate the 
integral over z, which yields the Sine function representa- 
tion of the propagator, 



GAh{x) 
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The massless limit of equation (||) is obvious, and setting 
m = makes no difference to its impressive convergence 
properties Q|. The difference between this form of Gah 
and the form found in Ref. Q is largely aesthetic: when 
m =/= 0, rescaling s to s/m makes all the terms inside the 
exponentials explicitly dimensionless. Lastly, it will be 
convenient to write G\h as 



Gjy h (a - b) 



(4tt) 2 



^2 p(fc)exp 



4c(k) 



where 
c(k) 
p(k) 



„kh 
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exp (kh — m 2 e kh ) 
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Feynman diagrams are integrals over products of prop- 
agators, and we derive the Sine function representation 
for a given diagram by inserting G\h for each line. All 
the spacetime integrals required by the internal vertices 
can be performed analytically, and we are left with an 
TV-dimensional infinite sum which approximates the orig- 
inal integral, where N is the number of internal lines in 
the diagram. In Ref. [| we give Sine function Feynman 
rules which describe the construction the Sine function 
representation for an arbitrary diagram, and these are 
not changed by the inclusion of the massless propagators. 

3 The Master Diagrams 

The 10 three-loop master diagrams, or vacuum tadpoles, 
are depicted in Figure In addition to their intrinsic 
beauty, their practical use is that three-loop diagrams 
dominated by a single, large mass and with small exter- 
nal momenta can be expressed in terms of the master 
diagrams. Four of these diagrams were given analytically 
before 1998 (q-fj, and Broadhurst found results for the re- 
maining six [0 . Applying the coordinate space Feynman 
rules to the overall "tetrahedron" topology of the master 
diagrams yield the following integral 

V, (A) = (4tt) 6 J d 4 x 1 d 4 x 2 d 4 x 3 G 1 (x 1 - x 2 ) 

G 2 (x 1 - x 3 )G 3 (xi - x 4 )G 4 (x 2 - x 3 ) 
G 5 (x 2 - x 4 )G 6 (x 3 - Xi). (10) 

The subscript j labels the specific combination of mas- 
sive and massless propagators, as enumerated in Figure [l|, 
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Figure 1: The three-loop master diagrams, where dashed 
lines are massless. 



while the prefactor reconciles our measure with Broad- 
hurst's. Equation (|l0| ) is a coordinate space integral, but 
this analysis works equally well in momentum space. The 
vertex labels are shown in Figure |[ and from here on we 
use Lorentz invariance to set a; 4 to zero. The six prop- 
agators are functions of the mass, rrij, carried by each 
line. For the master diagrams, the rrij are either or 1, 
but the Sine function representation yields results even 
if all the non-zero masses are different. The integral Vj 
is divergent, but it is regulated by employing the cut-off 
propagator, Ga(x). 




Figure 2: The generic topology of the master diagram, 
showing our choice of vertex labels. 

Applying the Sine function Feynman rules, we insert 
G\h for each of the propagators, and perform the Gaus- 
sian integrals. After some algebraic manipulation, we ob- 
tain the the Sine function representation of the general 
integral Vj , 



Vj(A,h) 



nLi c »(M 2 P"(M 

(a + 6c 6 (fc 6 ))2 ; 



(11) 



where 



a = (c 2 (fc 2 ) + c 5 (k 5 ))c 3 (k 3 )c4,(k 4 ) + 
C2{k 2 )c b {k 5 )(cz(k 3 ) + Ci{ki)) + 

Cl(fcl)(c 2 (fc 2 ) + C3(fc 3 ))(c4(fc4),+C5(fc 5 )) (12) 

b = (ca(fca) + c 4 (k 4 ))(c 3 (k 3 ) + c 5 (k 5 )) + 

ci(k 1 )(c 2 (k 2 ) + c 3 (k 3 ) + c 4 (fc 4 ) + cb(&b)).(13) 

The Cj and pt are defined by equations (@) and (^) and the 
summation over k is shorthand for infinite sums running 
between ±00 for each fc,-, i = 1 • • -6. 

To compute Vj from the Sine function representation 
we numerically evaluate the six-dimensional infinite sum 
in equation ([Ll]). In practice, Vj(A,h) is an approxima- 
tion to Vj(A), which is exact in the limit h — > 0]. By 
choosing h, we effectively fix the accuracy with which it 
is possible to evaluate Vj. Moreover, although the cutoff 
renders Vj finite, the result of evaluating equation jll| ) 
directly depends on A, and we must remove this regular- 
ization dependence before proceeding. 

Examining the six-dimensional sum, we see that if 
A — > and ki — k, then the general term approaches 
a constant, /i 2 /(256(47r) 6 ), as k becomes large and nega- 
tive. Consequently, the sum diverges in the limit A — > 00, 
and we have thus identified the origin of the divergent 
part of Vj in the Sine function representation. 

In Ref. @, we renormalized the Sine function repre- 
sentation of the sunset graph by expanding it about its 
external momentum and subtracting the p° and p 2 terms. 
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Table 1: We have computed the 9 non-trivial Fj and com- 
pared them to the exact values. These results were ob- 
tained with h — 0.7. 



This approach cannot be used with the master diagrams 
as they do not carry external momenta. However, the di- 
vergent part does not depend on any of the masses in the 
diagram. Consequently, the difference of any two Vj (A, h) 
is finite even when A — > oo, as the divergent parts in the 
sums for each of the Vj(A, h) will cancel term by term. In 
an approach that parallels Broadhurst's analytical treat- 
ment of the Vj, we write the differences between Vj and 
Vo as 

Fj(h)=J2} im ( 14 ) 

^ — * A^oo 
fci 

where tj is the general term from the six-dimensional sum 
that gives the Sine function representation of Vj ■ By form- 
ing the difference between V\ and Vj in this way, and then 
taking the limit A — > oo, we compute the finite Fj without 
ever needing to evaluate the Vj(A), which diverge in the 
limit A — > oo. By taking this approach we avoid the loss 
of precision that would result from attempting to evalu- 
ate Fj by computing Vj at large but finite A and then 
subtracting the result from V\. 

In Figure |l| we have used the symmetry of the dia- 
grams to ensure that G\, the propagator between x\ and 
X2, is always massive. Since all the other lines in V\ are 
massless and pk with m 2 > is always less than p^ with 
m 2 = 0, our choice of labeling ensures that t\ > tj for 
any {k\, ■ ■ ■ , kg}. Thus all the terms in Fj are positive, 
and a numerical evaluation of the multidimensional sum 
approaches the limiting value monotonically. 

We could be content with simply calculating the dif- 
ferences between diagrams and comparing our result to 
the analytic values. However, we can make use of a sin- 
gle analytical result, namely that the finite part of V\ is 
6((3) + 3£(4), to relate our Sine function representation 
calculation of the Fj with the renormalized Vj . 
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Figure 3: We plot 1 — Fio(h)/Fio, as a function of h. 

4 Numerical Evaluations 

We evaluate the Fj with a code that is similar to those 
used to obtain the results given in Ref. The code 

adaptively limits itself to the region of the six-dimensional 
lattice defined by the ki that make the strongest contri- 
bution to the overall sum. We accelerate the convergence 
of each of the six sums by using an Aitken S 2 extrapola- 
tion |j|,|lo||. This can significantly reduce the time needed 
to evaluate the sums, but the problem is still tractable 
without the extrapolation. While we have endeavored to 
ensure that the algorithm is implemented efficiently, there 
is no guarantee that it is the optimal method for evaluat- 
ing the Sine function representation of Fj . Consequently, 
the timing data we present (obtained using a 250MHz Ul- 
traSPARC II CPU) should be interpreted as a upper limit 
to the effort needed to evaluate these diagrams' Sine func- 
tion representations. 

When using equation (Q) to evaluate the Fj, we find 
that there are two distinct sources of error. The first is 
the intrinsic difference between Fj(h) and the exact value 
of Fj, and the second is the numerical error (if any) in 



. 1 Sa.mpln mHns for calculating Eg rart hp found at the URL 
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the evaluation of the infinite series. This has two possi- 
ble origins: either because the summation was truncated 
prematurely, or because of an accumulation of numerical 
error. In Table § we show the results for the nine non- 
trivial Fj, with h = 0.7. It is comparatively easy for us 
to obtain results accurate to within a few parts in 10 7 . 

The difference between the Sine function representa- 
tion and the exact value decreases quickly as h is reduced. 
In Figure ^, we plot the error in the value of Fio computed 
from the Sine function representation, as a function of 
h. The effort required to evaluate the sum increases at 
smaller h. First, the general term dies away more slowly 
when h is decreased, so we must add more terms to eval- 
uate the Sine function representation to a given level of 
accuracy. In addition, with smaller values of h, the maxi- 
mum size of the terms which can be ignored by the numer- 
ical algorithm also drops, which further increases the ef- 
fort needed to compute the multi-dimensional sums. With 
h ~ 1.0, the Sine function representation could be eval- 
uated in approximately three minutes on a workstation. 
The exact timings depend on the truncation procedure 
used in our code, which is somewhat heuristic and prob- 
ably too conservative. With h ~ 0.7, the evaluation time 
rises to around 20 minutes, and with h ~ 0.5, several 
hours are required. 

By paying close attention to the form of the sum, it 
is possible to improve on these speeds. We are particu- 
larly interested in computing Fq, the master diagram in 
which all lines are massive. A direct analytical treatment 
of this diagram yields a double integral, which must then 
be computed numerically in order to verify the results ob- 
tained with the PSLQ algorithm 0] . All the other results 
identified this way can be compared to (at worst) one 
dimensional integrals, which could be computed to very 
high precision. However, the result for Fq has not been 
checked beyond 15 decimal places. 

We can make two improvements to the efficiency of the 
Sine function representation of Fq. First, we can compute 
V5JV — Vq, instead of Fq, as the sums in the positive ki 
directions for i = 2, ■ • ■ , 5 converge much more rapidly 
than those in Fg, due to the exp (— e ' ) terms contributed 
by the extra massive propagators in Vgjv- Since Vsjv is 
known analytically, we can then deduce both V$ and F&. 
Second, we can speed the evaluation of the inner loop 
since we can write equation (|ll|) as 



where 



V 3 (A,h) 



E 



n n= i c n (k n )^p n {k n ) 
62 



(15) 



{a/b + c 6 (k 6 )y 



(16) 



and we recall that b does not depend on k§. Since f(a/b) 
depends on the single parameter a/b we can construct 
an interpolation table which allows us to evaluate / in 
far less time than it takes to compute the sum over k&- 
This effectively reduces the problem to a five-dimensional 
sum. The combination of these two methods improves the 
evaluation time by an order of magnitude or more, and it 
is straightforward to calculate Vq to within a few parts in 
10 13 . 

Unfortunately, going beyond 16 digits of precision re- 
quires the use of quadruple precision variables, and this 
dramatically lowers the computational efficiency of the 
code on present hardware. However, given a machine that 
performed 128-bit numerical arithmetic at similar speeds 
to 64-bit arithmetic, we would be able to evaluate Fq to 
a precision much better than 1 part in 10 -16 . 

5 Discussion 

The immediate purpose of this paper was to show that the 
Sine function representation reproduces the analytic re- 
sults for the three-loop master diagrams. Since our results 
are limited only by the finite accuracy of double-precision 
computer arithmetic, we have emphatically demonstrated 
the efficacy of the Sine function representation in this spe- 
cific case. We pursued this problem to test the usefulness 
of the Sine function representation as a tool for evaluating 
more general Feynman integrals. The Sine function rep- 
resentation does not rely on any special properties of the 
diagrams being evaluated, and we made no use of the ana- 
lytic knowledge gained from previous work on the master 
diagrams. Thus, these results support our claim that the 
Sine function representation may be a useful approach to 
evaluating general higher order diagrams Q . 

There is considerable theoretical interest in perform- 
ing accurate calculations of specific diagrams to assist 
the investigation of non-trivial relationships between dia- 
grams. Broadhurst's evaluation of the three-loop master 
diagrams is a prime example of this type of work, which 
has the ability to illuminate the deep structure of pertur- 
bative quantum field theory. The Sine function represen- 
tation has the potential to facilitate this approach, since it 
is fast, accurate, and does not require a partial analytical 
evaluation of the diagram. 

When analyzing experimental data, we are unlikely to 
need the high levels of accuracy achieved in this paper. 



5 



However, if we only want a few significant figures, we 
can increase h, which allows us to evaluate diagrams very 
rapidly. The Sine function representation's convergence 
properties, lack of analytic overhead, and applicability 
to arbitrary topologies suggests that it may be a useful 
tool for automatically calculating large sets of diagrams. 
However, before this possibility can be explored in de- 
tail, we must generalize the Sine function representation 
to fermionic and vector fields, and describe a renormal- 
ization procedure that can be applied to the Sine function 
representation of an arbitrary diagram. This work is cur- 
rently in progress. 

Finally, while Feynman integrals arise in perturbative 
quantum field theory, the Sine function representation 
is an outgrowth of a new approach to "exact" numer- 
ical quantum field theory, the Source Galerkin method 
^l|-|l6|. To apply the source Galerkin method we must 
carry out integrals over products of terms proportional to 
the free propagator, which resemble Feynman integrals. 
Consequently, the numerical codes and analytical insight 
needed to apply the Sine function representation to Feyn- 
man integrals also improve the efficiency of the Source 
Galerkin method. 
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